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Abstract 

We present an improved version of our program package oneloop which - written 
as a package for MAPLE jl], || - solves one-loop Feynman integrals. The package is 
calculating one-, two- and three-point functions both algebraically and numerically to 
any tensor rank. In addition to the original version oneloop 2.0 also calculates infrared 
divergent integrals. Higher powers of propagator terms and the 0(e) parts relevant 
for two-loop calculations are now supported. 
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NEW VERSION SUMMARY 



Title of new version: oneloop.ma, version 2.0 



Reference to original program: oneloop.ma 

Reference m CPC: Comput. Phys. Commun. 85 (1995) 153 

Does the new version supersede the old program? Yes 



Licensing provisions: None 



Computers: Any platform which has MAPLE V, Release 3 



No. of processors used: 1 



Operating system under which the new version has been tested: (i) Unix (Linux 2.0) (ii) VMS 6.2 
(iii) MS-DOS 6.0 



Programming language: MAPLE V, Release 3 

Subprograms used: cfcn.ma, r.ma, simple. ma, pv.ma, mess. ma 



Memory required to execute with typical data: Less than 8 MB RAM 

Typical running time: The algebraic or numeric evaluation of a three-point function up to the 
rank-3-tensor on a DEC Alpha 8400 needs less than 2s 

No. of lines in distributed program including subprograms: 4019 

No. of bytes in distributed program including test data: 168 189 

Distribution format: binary (zipped ASCII files) 

Keywords: Electroweak theory, Feynman diagrams, one-loop corrections, renormalization 

Nature of physical problem: The theoretical determination of cross sections in particle processes 
requires the calculation of radiative corrections. The most important contribution comes from the 
level of one-loop Feynman diagrams which arise from the model under consideration, usually the 
standard model of elementary particles 

Method of solution: This package is designed to evaluate automatically one- loop integrals. It 
is making use of the properties of 1Z functions, a class of special functions which simplifies the 
evaluation of Feynman integrals 



LONG WRITE-UP 



1 Introduction 



In particle physics the calculation of one- and higher loop corrections to particle processes 
is mandatory to keep track with the increasing accuracy of particle colliders. For this 



reason the necessary calculations were automated by virtue of computer programs during 
recent years. In this context the program oneloop was developed at the one-loop level pi. 
It provides the necessary integrals which are needed for Feynman diagrams with at most 
three external legs. 

Several improvements of the original version now merged in a new version 2.0. It 
contains the following new features: 

• Arbitrary powers of the propagator terms in the denominator are allowed (cf. sect. ^). 
The original version was restricted to propagators of power one. 

• Infrared divergences can be calculated in dimensional regularization. They will ap- 
pear as poles in the dimensional regulator e = (4 — D)/2 like in the ultraviolet case 
(cf. sect. |). 

• The contributions of 0{e) can be calculated as well. Although these terms are not 
contributing in pure one-loop calculations they become relevant in two-loop appli- 
cations (cf. sect. H). 

• Several kinematical points in which the original version ran into numerical instability 
or terminated in a division by zero are now solved in a different, numerically stable 
way (cf. sect. ||). 

• Additional functions OneLoopTensnPt accept squared momenta as input and return 
full tensors (cf. sect. |6|). 

• The usage of a library of one- loop integrals is generalized (cf. sect. ^). 

• Some other minor changes have to be reported (cf. sect. ^): 

In the case of the three-point function the convention of the output changed slightly. 
It has now the same structure as for the other functions. 

Some new functions which abbreviate the output are introduced in the case where 
the functions are calculated analytically. 

This note describes all changes in more detail. 

2 Arbitrary powers of propagators 

In addition to the original functions 

• OneLooplPt(p, m) 

• OneLoop2Pt(j>o,£>i, 9, m i, m 2) 

• 0neLoop3Pt (p , p 1 , p 2 , qi , q 2 o , 921 , mi , rn 2 , m 3 ) 

- which need the tensor degree, the external momenta q n and the masses m n as input - 
there exist now the following generalizations 

• OneLooplPt(p, m, t) 

• 0neLoop2Pt (p , Pi , q, mi , m 2 , t\ , t 2 ) 

• 0neLoop3Pt(p ,Pi,P2, 91,920, 921, mi, m2,m 3 ,ti,t 2 ,*3)- 



These functions additionally expect the powers t n of the propagator terms in the denom- 
inator. In detail they directly correspond to the following integrals: 



• One-point function: 

OneLooplPt(p,m,t) = A^\m) = J d D l „ 2 _ ( ^" (1) 



• Two-point function: 




OneLoop2Pt(po,pi, q, m 1 ,m 2 ,t 1 ,t 2 ) = B {poPl){tlt2) (q, m 1 ,m 2 ) 



[(/ + q) 2 -m\ + iqY* [I 2 - m\ + ig] t2 



Abbreviations: 



l-q 



2 ' 



i ± = q - p 



Three-point function: 




(2) 
(3) 



0neLoop3Pt (p , pi , p 2 , qi , q 20 , qi\ , mi , m 2 , m 3 , h , t 2 , t 3 ) 
= C( poPlP ^' t2 ' t3 \ qi ,q 2 , mi ,m 2 ,m 3 ) 

(i \\) po (h\\r v±r 



i 



d D i 



(4) 



[(/ + qi ) 2 -m\ + ig]^ [(I + q 2 ) 2 -m\ + ig] 1 * [I 2 -mj + ig]* 



Abbreviations: 



l o\\ = —j= ; = F^T ; 9 2 = 92 — 9i (5) 

h = Jilt - In - I 2 ; 920 = ^-3 ; g 2 i = \l gfo - 5§ 



Our notation for the OneLoopnPt functions distinguishes between parallel and orthogonal 
space which is reflected by the definitions of (g) and (g) for the momentum components 
U|, |5|. Please keep in mind the fact that in this notation the indices Po,pi, ■ ■ ■ represent the 
powers of the different components of the loop momentum I which should not be mixed 
with Lorentz indices. In our notation hi, l n, l^n describe the components of I which are 
parallel to the external momenta q,qi,q2, whereas l± represents the orthogonal comple- 
ment. 

If the powers t n of the propagator terms in the denominator are omitted the program 
assumes the default value 1 for each t n . 



3 Two-loop relevant parts of one-loop diagrams 

The results of the OneLoopnPt functions are Laurent expansions in terms of the ultraviolet 
regulator e. Therefore the output of the OneLoopnPt procedures consists of a list where 
the significant coefficients (C(e _1 ), O(e )) of this expansion are given. 

In two-loop calculations one-loop diagrams get multiplied with divergent Z-factors. 
Two-loop integrals can factorize into a product of one-loop integrals. This is the reason 
why one- loop contributions of 0{e l ) are also of interest |], f?J. For that purpose this 
coefficient is calculated if the qualifier more is used in the function call: 



input 


output 


OneLooplPt(p, m, t) 


[e~ 1 -term, e^-term] 


OneLooplPt(p, m, t,more) 


[e~ L -term, e^-term, e 1 -term] 


0neLoop2Pt (p , pi , q, mi , m.2, t\ , t 2 ) 


[e~ 1 -term, e^-term] 


OneLoop2Pt(po>Pi 5 q, mi, m2,ti,t 2 ,moTe) 


\e~ x -term, e^-term, e 1 -term] 


0neLoop3Pt (p , Pi , p 2 , qi , q20, 921 , m\ , m 2 , m 3 , t\ , t 2 , h) 


[e~ l -term, -term] [ 


0neLoop3Pt (p , pi , P2 , 91 , 920, 921 , 

m 1 ,m 2 ,m 3 , t\, t 2 , i 3 ,more) 


[e~ 1 -term, e°-term, e 1 -term] 



In the output "e _1 -term" is meant to represent the divergent part (the coefficient of 1/e) 
whereas "e°-term" is describing the finite part and "e 1 -term" the 0(e) contribution. 

As long as the arguments of the OneLoop functions are symbols the output is given 
algebraically, whereas numbers inserted for the arguments imply a numerical result. 



4 Infrared divergences 

Infrared divergences are now regulated as well. This kind of divergence may occur for 
instance for the three-point function if the momenta are on-shell and one mass is set to 0. 

lr The output of 0neLoop3Pt changed slightly compared to the original version (cf. sect. H). 



This effect also appears for a two-point function where a propagator term with vanishing 
mass is squared. Since the collinear divergent case is still excluded the divergence will 
always occur as a pole term of 0(e^ 1 ). It means that the output involves a non- vanishing 
e~ 1 -term in the notation of sect. ||. There is no explicit distinction made between UV and 
IR divergences, there is no infrared dimension parameter sir,. Both kinds of divergences 
are described by e = (4 — D)/2. The result for the collinear case is well-known It is 
the only case which has a pole of 0(e~ 2 ). We excluded this simple case, but the user can 
add it easily himself. 



5 Numerical stability 

The accuracy and stability of numerical results may suffer from cancellations of large, 
approximately equal dilogarithms. Since MAPLE supports calculations of arbitrary pre- 
cision, it is possible to increase the number of digits to improve accuracy. Usually we 
calculated with 20 or 40 digits. 

Nevertheless, at some kinematical points the general procedure breaks down since 
divisions by or similar errors are encountered. It turns out that these points always 
belong to a kinematical arrangement - for instance equal or vanishing masses or momenta 
- which is simpler than the general case. For this reason, even if there were no numerical 
problems - in the sense of optimizing the code - it is advisable to solve the simpler 
kinematical situation in a faster than the general way. 

This is done completely automatically for the following kinematical points: 



0neLoop2Pt: 

9 = 
mi = 

7772 = 

(7 = , 777l = m2 

q 2 - m\ + ml = 



q 2 + ml — ml = 
q + 777i + m 2 = 
q — 777i + m 2 = 

q + 777l — TTi2 = 

q — mi — m% = 



0neLoop3Pt : 

qi = 
921 =0 

7771 = 

777 2 = 

777 3 = 

920 - 91 = 921 
920 - 91 = -921 
920 = 921 
920 = -921 

qi = 0, ?77l = 777 3 

920 = 0, 921 = 

921 = 0, 9i = 920 

921 = 0, (q% - 777^)91 - (ql - ml)q 20 

- mf(9 2 - 91) = 
920 = 0, 92i = 0, 777 2 = m 3 



921 = 


= 0, 91 


= 920, mi = 


777 2 








9i = 


0, 777i 


= m 3 , \/q\ Q 


-921 


+ 777l - 


\- 777 2 


= 


9i = 


0, 777-1 


= ma, y/q$ 


-921 


— 777l ~ 


1- 777 2 


= 


9i = 


0, 777-1 


= m 3 , a/ 920 


-921 


+ 7771 " 


- 7772 


= 


9i = 


0, 77ll 


= ma, y/ql 


-921 


— 777l " 


- 7772 


= 



920 = 


= 0, 921 = 0, 777 2 = 


7773, 91 + 7771 + 777-3 : 


= 


920 = 


= 0, 921 = 0, 7772 = 


m 3 , 91 - 


- 777-1 ~ 


h 777-3 : 


= 


920 = 


= 0, 921 = 0, 777 2 = 


m 3 , 91 4 


- 777-1 " 


- m 3 = 


= 


920 = 


= 0, 921 = 0, 7772 = 


m 3 , 91 - 


- 777-1 " 


" "7-3 = 


= 


921 = 


= 0, 91 = 920, mi = 


= m-2, 91 


-1- 7771 


+ m 3 


= 


921 = 


= 0, 91 = 920, mi = 


= m-2, 91 


— 7771 


+ m 3 


= 


921 = 


= 0, 91 = 920, mi = 


= m 2 , 91 


f 777l 


- 7773 


= 


921 = 


= 0, 91 = 920, mi = 


: m 2 , 91 


— 777l 


- m 3 


= 


9?- 


777i + m\ = 










q'l + 


ml — ml = 










9?- 


m| + m§ = 










ql + 


777-I — 7773 = 










ql + 


777^ — 7772 = 










q\- 


777i + 7772 = 










920 - 


- 7772 + m 3 = 










<&o- 


-ml — ml = 












- 921 - m l + m 3 = 











9lo - 


- 92i + m| - m 2 = 











(920 


— 9l) 2 + 777i — 7772 


= 








(920 


— 9l) 2 — 777i + 7772 


= 








920 - 


- 921 - mf + ml = 












120 ~ ill + m l -m\ = Q feo - <7i) 2 - ?Ii = m v ?Io - ill = TO 3> m 2 = 

920 - 921 = m 2> Ql = m l, m 3 = feo - qi) 2 - ill = m 2: 1l = m l> m l = 

Several of these points correspond to different thresholds of the integrals. There are 
only two remaining sources of difficulties: 

• if an integral is first solved in a general manner and then - in the result - special 
values are substituted which correspond to one of the aforementioned kinematical 
points. Solution: the values have to be assigned from the beginning of the calculation. 

• if the kinematical arrangement is not in, but very close to (less than roughly 10~ 10 
compared to the relevant scale) one of the above listed kinematical points. Solution: 
the value of Digits has to be increased. 



6 Tensors 

The functions OneLoopTensnPt expect squared momenta - q 2 for two-point functions, 
q\,q\ and (q 2 — Qi) 2 for three-point functions - as arguments and return a full tensor 
as output - which is completely equivalent to the notation of (||) and (^) but perhaps 
the more familiar representation ||. They replace the similar PassVelt functions of the 
original version which now are no longer needed. 

The following types of arguments - here only demonstrated for the two-point case - 
are allowed for all OneLoopTensnPt functions: 

• 0neLoopTens2Pt(i): 

The rank i tensor decomposition of the two-point function is given. The procedure 
returns a list consisting in different coefficients, which are expressed in terms of 
the OneLoopnPt tensor integrals, and the defining equation for the coefficients, for 
instance in the case i = 2: 

0neLoop2Pt(0, 2) DneLoop2Pt(0, 2) 

^21 — 7. 7, ) ^20 — o7 — o 1 o \ — 

3 — 2e q l (—3 + 2s) 

0neLoop2Pt(2,0) 

H ^2 ' ^20^1^2 + ^2i57tiM2 

0neLoopTens2Pt(z, full): 

The function inserts the results of the OneLoopnPt tensor integrals explicitly. Here 
again the case i = 2: 



C 2 i = [e 1 -term, e° -term] , C 2 o = [e 1 -term, e Q -term] , C 2 oq^ 1 q fl2 + C 2 ig fil 



l>2 



• 0neLoopTens2Pt(i, q 2 , mi, 7712, ^1,^2): 



The function returns the same as 0neLoopTens2Pt(i, full), but expressed in the 
user defined terms for q 2 ,m\ and m^. Of course numerical values are also allowed. 
t\ and t<i are optional. 

0neLoopTens2Pt(i, q 2 , m\,m2, t\, £2, more): 

The function returns also the 0(e) contribution. Again, t\ and t 2 are optional: 

C21 = [e -1 -term, -term, e 1 -term] , C20 = [e~ 1 -term, e° -term, e 1 -term] , 
C20^i?M2 + C2i9niim] 



7 Library 



The usage of a library of integrals was restricted to the PassVelt procedures in the original 
version. It is now applicable for all OneLoop and OneLoopTens procedures, so that in any 
case the procedures speed up. The package now includes the functions 

• 0neLoopLib2Pt(i, j) 

• 0neLoopLib3Pt(i, j) 

which generate a library of all 0neLoop2Pt and 0neLoop3Pt functions respectively up to 
the tensor rank i. All powers of the denominators from 1 to j are considered. The second 
parameter j is optional. If it is omitted the procedures assume the value 1 for j, so that 
no powers of denominators higher than one are calculated. 

If this library shall be written in any other than the current directory one has to assign 
the variable LibPath which substitutes the variable tensorpath in the original version 

LibPath : =(path) ; 

Now the library will be written to (path). The procedures will only look for the library if 
LibPath is assigned. If the procedures search for an integral which is not contained in the 
library, an error message is returned. 

In practice it is necessary to store not only the general mass case of each integral. 
Several special cases mentioned in sect. |5| are also needed. The necessary integrals are 
automatically created by the OneLoopLib routines. Each integral corresponds to one file. 



8 Minor changes 

We changed the convention for the three-point function compared with the original version: 
The vanishing e~ 2 -term which was kept for testing reasons is dropped now. The list of 
abbreviations which came with the three-point function also disappeared, because it is no 
longer needed. 

Instead the output now is written in the algebraic case by using several new abbrevi- 
ations: 



R2exl(x, y) 



1 



x 
V 



In 1 



1 



In 1 + Jl 



+ in 



(6) 



ln(— x) 



ITT 



R2ex2(x,y) 




1 - £ 

y 



1 + ,/1-f 



+ 1 + 



+ 1 



+ 1 





X 






y 




X 








y 




X 








V 



In 1 



1 2 



1 + ,/1-f 



1 - 2 

y 



In ( 1 + ,/!-- 

y 



(7) 



+ — (lny) 2 + 2 (lnx) 2 — 2 lnx lny 



(In y — 2 In x) 



1 + W1 



In 1 



1 



1 



In 1 + 



x 



-in ^LJt [2 ln2 + ln(x-y)] 



-y 



R3ex2(x,y,z) = 2 In ( 1 - r/(x, 2) + 2 In ( 1 - ~j 7?(y, z) 
-1-2 Li 2 ( 1 - -) + 2Li 2 f 1 - ^ + 2 (lnz) 2 



(8) 



R3ex3(x, y, z) 



4S 12 I 1 - -) - 4Li 3 { 1 - -1 +4Si 2 ( 1 - ^ - 4Li 3 f 1 - ^ 



4 77 x, 



1 



1 



Li 2 ( - ) + In ( 1 

z, 



4^2,,- Li 2 i +ln 1-i In 



In 



x\ 1 
z) 2 



In 1 



In ( 1 - - 

z 



1 2 



— Alnz 

- 4 In 2 
+2 Li 3 



.7; 



Li 2 1 - - + Li 2 1 - - 



/? ■ x ' ~) ln v 1 ~ / + v \ y " z) ln v 1 ~ ~ 



+ 



In 



z-y 

z — X 

z-y 

z — X 
' l In 



n 2 



( (z — y)x\ T / x 

-2Li 3 ) ^- +2Li 3 - 

\(z-x)yj \y 

( y — x z \ ( 1 

77 , ) +2rj[y-x,- 

\ z z — X \ z 



+2 (lnx — lny) 



z-y 



z — X 



Li 2 



ln 



[z - y)x 



z — X 



z — y\ l z — y z 
rj 



z z — X 



+2 



(z - x)y 

y?-'< 

lnx — lny In ( 

2 \(z-x)y 



Li 2 



(i) 



In (lnx-lny) 2 -2C(3)-^(lnz) 3 



(9) 



- x)y 



y z — x 



These functions are related to the corresponding 1Z functions^] |T^, 11]. They represent the 
coefficients of the Taylor expansion in e making use of (cf. J12|]) 



C(n) 



00 -1 

y — 



fc=l 



2 The function R3ex3(a;, z) assumes that a; and y have an imaginary part of different sign which is 
always the case. 



Li2(2 ) _ -j^^ds 



Li 3 (.) = j^ds 



1 / ln 2 (l- S ) 

s 12 (z) = _y — — ^ 



7](a, b) = 2iri [0(-Im a)0(-lm 6)0(Im (ab)) - 0(lm a)0(lm 6)0(-Im (ab))\ . 
The function 
evalRex(( expression)) ; 

substitutes the Rnexm functions in (expression) by the corresponding Logarithms and 
Dilogarithms. 

9 Run time 

Due to several optimizations in program code the procedures became much faster compared 
to the original version. Run times typical for the new version are displayed in the following 
table. We used three different systems: (i) a DEC Alpha 8400 workstation (ii) a VAX 
4000/90 workstation (hi) a PC with i486 chip, 16 MB RAM and 66 MHz frequency All 
systems are working with MAPLE V. Without using the library we found the following 
run times: 





DEC Alpha 8400 


VAX 4000/90 


PC/DOS 


tensor degree 


numerical 


algebraic 


numerical 


algebraic 


numerical 


algebraic 


two-point functions 





0.02 s 


0.04 s 


0.11 s 


0.22 s 


0.50 s 


0.88 s 


1 


0.04 s 


0.21 s 


0.49 s 


1.53 s 


1.64 s 


2.70 s 


2 


0.07 s 


0.28 s 


0.50 s 


1.97 s 


1.79 s 


3.90 s 


3 


0.08 s 


0.41 s 


0.57 s 


2.99 s 


2.14 s 


8.51 s 


three-point functions 





2.54 s 


1.83 s 


19.16 s 


13.82 s 


92.5 s 


25.87 s 


1 


2.57 s 


2.23 s 


20.68 s 


16.66 s 


117.8 s 


26.36 s 


2 


3.01 s 


3.62 s 


20.71 s 


26.31 s 


120.9 s 


38.89 s 


3 


3.78 s 


4.16 s 


28.87 s 


40.30 s 


125.7 s 


117.7 s 



It should be emphasized that these are the times which are necessary to generate the 
functions once and forever using the OneLoopLibnPt routines. If they once are stored in 
the LibPath directory they may be read in quickly. With assigned LibPath we get: 





DEC Alpha 8400 


VAX 4000/90 


tensor degree 


numerical 


algebraic 


numerical 


algebraic 


two-point functions 





0.34 s 


0.03 s 


2.40 s 


0.09 s 


1 


0.39 s 


0.05 s 


2.95 s 


0.15 s 


2 


0.43 s 


0.06 s 


3.22 s 


0.19 s 


3 


0.47 s 


0.06 s 


3.43 s 


0.21 s 


three-point functions 





1.30 s 


0.09 s 


9.27 s 


0.45 s 


1 


1.32 s 


0.10 s 


9.97 s 


0.48 s 


2 


1.36 s 


0.12 s 


10.30 s 


0.56 s 


3 


1.88 s 


0.45 s 


12.08 s 


0.65 s 



Appendix: Installation and program call 

The distribution of the oneloop package is accessible at 
http : / / wwwthep . physik . uni-mainz . de/~xloops/ 

To install the package, just copy all files with the extension .ma in one directory. After 
having started the MAPLE session, you can read the oneloop package - just type the 
following lines: 

LoopPath : =(path) ; 
read' (path) one loop .ma' ; 

(path) describes the directory where the .ma files are located. If these files are in the current 
directory then you can skip the assignment of LoopPath and only type 

read' oneloop .ma' ; 
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